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Abstract 

Stable assemblages of localized vortices exist which have particle-like 
properties, such as mass, and which can interact with one another when 
they closely approach. In this article I calculate the mass of these localized 
states and numerically investigate some aspects of their interactions. 

(N 

^ ! Introduction 

qv I The system of N point vortices on a sphere is a Hamiltonian system with symme- 

^\ • try 5*0(3). Steadily rotating arrangements of vortices are the relative equilibria 

of these systems, and the nonlinear stability of these relative equilibria may 
be analyzed as a matter of routine, at least in principle. If stable then such 
an arrangement of vortices will, under sufficiently small perturbation, approx- 
imately maintain its "shape" , so that its subsequent motion will be nearly its 
initial condition up to the action of the symmetry group 5*0(3). However, if 
the momentum of the relative equilibrium (i.e. its center of vorticity) is zero 
then the position of such a relative equilibrium on the sphere is not stable under 
perturbation to nonzero momentum, as has been established in Patrick [1992, 






X 



1995, 1998[ 



In this article I show that there are arbitrarily localized relative equilibria of 
N vortices with zero momentum, and that these relative equilibria are formally 
stable for N — 3 and N = 4. Under small perturbation to nonzero momentum, 
they move on the sphere, to first order in the momentum perturbation, as a free 
particles move under the influence of a magnetic monopolc. The effect of the 
monopole diminishes as the relative equilibria become more localized. Since a 
given momentum perturbation implies a particular velocity of motion of each 
relative equilibrium as a whole, they have mass, and using a perturbation theory 



of Patrick [1995|, I find a formula for that mass. 



For N — 4, when two such relative equilibria are placed on the sphere, 
they will move independently until they approach within a distance comparable 
to their size, whereupon an interaction occurs. For example, one might collide 
with its opposite (obtained by reversing the first's component vortex strengths), 
the result being "radiation" in the form of dipole pairs, and generally opposites 
repel and likes attract. Here in this article I numerically determine some general 
features of these interactions, which are the subject of ongoing investigations. 

Over small regions of the sphere, the point vortex system on the sphere 
becomes the point vortex system on the plane, and since the above relative 
equilibria may be arbitrarily localized, one expects such in the planar case as 
well. Towards the end of this article I show this in indeed the case. The planar 
system of N point vortices will be a worthy motivating example for extensions 
of the stability theory of Hamiltonian systems with symmetry to the case of 
non-equivariant momentum mappings and/or noncompact groups. 

1 Context and Notations 

I briefly summarize here the basic elements of the system of N point vortices 



on the sphere and on the plane. For more details, see Kidambi and New 



bon [199J], Pckarsky and Marsden [1998], and the many references therein. 

The system of N point vortices of strengths r, S R on the sphere S 2 of radius 
R is the Hamiltonian system with phase space P sp = {S 2 ) N and Hamiltonian 

if sp (xi,-- • ,x N ) = y^ r„r m in/ m „ 2 , (i) 

7n<n 

where x n represents the location of the n th vortex and l 2 an = 2(R 2 — x m ■ x n ) is 
the chord length between x n and x m . The symplectic structure o/ p on P sp is 
a direct sum of the symplectic structures on S 2 weighted by vortex strengths; 
specifically 

N 

<^®lf^ ( 2 ) 



i R 

n—l 



where 



uj S 2(x)(v,w) = -—^x- (v xw). (3) 

Most significant from a geometric mechanics point of view is that this system ad- 
mits the symmetry of diagonal multiplication of 5*0(3) on P sp with momentum 
map 

1 N 

n—l 



These notations are as in Pckarsky and Marsden [1998]. 



In the system of N vortices in the plane the n th vortex has location z n = 
%n + i]Jn € -P pl = (C 2 )^, and the Hamiltonian and symplectic form are 



IP 1 



1 N 

-— ^2 r„r m in|z„ - z m \ 2 , w pl E0r n w o 



(5) 



where generally I will use the notation uj (a,b) = — Im(a&) for complex num- 
bers a,b G C This system admits the symmetry group SE(2) = {(e ie ,a)} of 
Euclidean symmetries acting diagonally on each factor C of P pl by (e l9 , a) ■ z = 
e l6 z + a, and a momentum mapping is 



N 



J pl 



E r ^ 



\W? 



(6) 



where se(2)* is identified with se(2) = {(0, a)} = R 3 by the standard inner 
product of R 3 . These notations are as in |Lewis and Ratiu 1996]. 

When attention is confined to a sufficiently small region of the sphere, the 
spherical system becomes the planar one. For example, using on each factor the 
map 



4>{z) = (z, ^r 2 - W) 

to pull back the sphere to the plane one obtains 



* ,sp _ 
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.P 1 



tjjsp 



R 2 



-IP 1 



(7) 



(8) 



Thus the pulled back equations of motion are approximately the same as the 
equations of motion of the planar system, since the factors of —1/R 2 here mul- 
tiply both the symplectic form and the Hamiltonian. 



As noted in Adams and Rath [198S], the momentum mapping of the planar 
system is equivariant if and only if the total vortex strength £ n T n vanishes. 
Indeed, the adjoint and coadjoint actions of SE(2) are 



Ad 



(e««,a)(y>a) 



/, e a — iffa), 



(9) 



CoAd( e ie !o) Gu,z/) = {ii-uj {e ie v,a),e w v) 1 



and the deviation of the momentum map J pl from equivariance is the cocycle 
a : SE{2) -> se(2)* defined by a = J pl ((e ie ,a) • p) - Coad (e * V) J pl (p) (the 
evaluation point peFis irrelevant). For the case at hand the cocycle is 



o-(^,a) = -(£ n r n ) 



\\a\ 2 
ia 



and the derivative of this cocycle at the identity is the skew-symmetric two form 
S : se(2) 2 -> R given by 



Z (01, ax) ,(62,0,2)) = (J2n T n)^o{ai,a 2 )- 



(10) 



Generalities on non-equivariant momentum mappings may be found in Abraham 



and Marsden [1978|; a prime fact is the momentum commutation identity 

{J £ ,J„} = J K) „ ] -E(&'7). (11) 

Both the planar and spherical systems have simple closed form solutions in 
the case N — 2: for the sphere any two vortices evolve as the action of the one 
parameter group with generator J sp /2ttIi2 while for the plane they evolve as the 
action of the one parameter group with generator 6 — (Ti+r2)/27r|zi — Z2I 2 , a = 
— J p1 /2tt\zi — Z2\ 2 ■ Numerically the action of these one parameter subgroups is 
easily computed, and since the Hamiltonians H sp and IP 1 are sums of pairwise 
interactions, the full system of N vortices may be numerically integrated in a 
symplectic, symmetry preserving and momentum preserving way using splitting 
methods, as in |Channcll and Ncri [ 1993 |. 



2 The Relative Equilibria 

To find the relative equilibria for the system of N vortices in the sphere one 
seeks £ e € so(3) ^ R 3 and p e e P sp such that dH sp {p e ) - dJ ie (p e ) = 0. Equiva- 
lently, using the obvious extensions of H sp and J sp to (R 3 )^, one can solve the 
equations 

W m H sp - V m J| c P = AV (\\x m \ 2 ) = A m x m , to = 1, • • • ,7V, 

where the A m £ K are Lagrange multipliers. Inserting ([j]) and (||) yields the TV 
equations 

P O77- f?2 

2J jp- = \ m x m + 27ri?£ e , A m = — - — A m , m = 1, •••,iV. (12) 

General analytic solutions are not to be expected to these nontrivial nonlinear 
equations. 

I seek a manifold of solutions to (|l2| ) contained in the zero level set of J sp , 
which in some limit is confined to arbitrarily small regions of phase space, and 
which is formally stable near that limit. I try a regular polygonal configuration 
where the first N — 1 vortices of equal strengths surround the (possibly different 
strength) iV th central vortex, N > 3. I will denote the strength of the central 
vortex by I\ For convenience I locate the central A^ th vortex at Rk where 
k = (0,0, 1), and the first vortex at (i?sina,0, i?cosa), so that a is the opening 
angle the angle at the center of the sphere between any outer vortex and the 
central vortex. The momentum is zero if and only if 

ri(xi +X2 H \-xn-i) +Tx N = (N - l)rii?cosak + ri?k = 

which is equivalent to 

r 

ri= ~(7V-l)cosa' (13) 



while equations (U2) reduce to Ai = • • • = \n-i and the two equations 



TiX + -g-k = Ai»i + 2TrRCe, (14) 



"1JV 
(N - l)IYRcOSa 

'ijv 



k = A A ri?k + 27ri?^ e , (15) 



where £ e depends on T, JV, a and -R, and where 



X 2 , X 3 Xn-1 



Setting 



fe=i 
so that 



^o = ? + ? + - + P- ( 16 ) 

'21 l 31 'JV-1,1 



W-l JV-1 

^ N = 2-^i e -27rik/N _ l ' Un = 2^1 | e 27rifc/Af _ ]|2 

fc— 1 ' 



M^ + vn = y^ 



v ' e 2TTik/N 

fe= 



_/ | e 27rifc/JV 1 i 2 ' 

one sees by symmetry that /ijv and v^ are real, and also that 

»n = -\{N-1), v N =^{N 2 -l). 
Then scaling /i/v-i and vjv-i gives 

Xo . 1= M*-i+»*-i Xo . j = 0; Xo . k= ^zii^. (17) 

it sin a i? sin a 

Equation (nq) implies £ e is along k and in (|l4| ) it is clear that the j component 
of both sides is zero, so writing (H4) and (|15[) in components gives three linear 
equations in the three unknowns Ai, \n, and £ e • k. These equations can be 
routinely solved to obtain 



*!=- „ * a . v ^o,r\/' . as) 



(JV-2)(JV-6) 
i? 2 cos a sin 2 a 12(iV - 1) 
T / iV 



Xn = ~ ~^ ■ 2 - ( COSQ! + ^77? — Tt)' ( 19 ) 



£ e -k = 5— ( hcosa). (20) 

5 47ri? 2 sin 2 a V^- 1 / 

In the limit a — > 0, from (Oh, the total vortex strength 

E„r w = r + (iv-i)r 1 = rfi--M (21) 



vanishes while the regular polygon of outer vortices collapses upon the central 
vortex. Although these relative equilibria depend on a, T and N, I will denote 
them below simply by p e . 

Formal stability means dehniteness of the Hessian of the reduced Hamil- 
tonian at the relative equilibrium, and, by the energy-momentum, method as 
in Marsden and Ratri| [1994], this is equivalent to dehniteness of the function 
H sp — J^ e on a subspace tangent to the momentum level set and complimentary 
to the tangent space to the group orbit. For N — 3 the reduced spaces are 
points and formal stability is immediate. I have determined the stability p e for 
N = 4, 5, 6 with the aid of a symbolic manipulator. For N = 4 the relative 
equilibria are formally stable without conditions, and in particular formally sta- 
ble arbitrarily near the limit a — > 0. For N = 5, but numerically now, one 
has formal stability if and only if a > 1.951 « 111.8° and for N = 6 if and 
only if a > 2.245 « 128.6°. Presumably this pattern continues and the relative 
equilibria are formally stable for N > 6 if and only if a is sufficiently near it. 

As already mentioned, since the momenta of such relative equilibria are 



zero, by Patrick [1992], formal stability implies stability only modulo 50(3); 
dynamically, these relative equilibria, when perturbed to nonzero momentum, 
approximately maintain their shape, which oscillates on a fast time scale, while 
they move around the sphere on a slow time scale. The left of Figure (Qj) 
is the result of a simulation and shows a typical motion of one such relative 
equilibrium. For short, I will call a point of phase space that results from 
perturbing one of the relative equilibria p e a "preq" . 

The motion of preq may be approximated as a direct application of the per- 



turbation theory of Patrick [1995|. In summary, the symplectic 50(3) symmetry 
implies that the linearization of the relative equilibrium has double eigenvalues 
and hence can be expected to have a nonzero nilpotent part, say N a . It can be 
shown that the image of N a is contained in the tangent space to the group orbit 
and the tangent space to the momentum level set is contained in the kernel of 
N a . Consequently there is a unique (it can be shown to be symmetric) bilinear 
form (here with the same name) N a on so(3)* such that the following diagram 
commutes 



T Pe P sp 



N a 



T Pc P sp 



dJ(pe) 



so(3)* 



N a 



?^«Pe 



50(3) 



I have calculated the N a in the case N = 3 and N = 4; they are the diagonal 
matrices with entries Jati, Jn2, Jn3 where 



J31 = J: 



1 



32 



87ri? 2 cos a ' 



1 



'33 



8ttR 2 sin a 



(2 cos 3 a + 3 cos 2 a + 2 cos a — 1 ) , 



(22) 
(23) 



and 

1 3 cos 2 a + cos a + 2 

41 = 42 = ~ 2irR 2 (1 - cos a)(9 cos 2 a + 4 cos a + 3) ' ( } 

J43 = j — (3 cos 3 a + 4 cos 2 a + 3 cos a — 2) . (25) 

YlirR 1 sin a 

The perturbation theory asserts that to first order the motion of the preq is 
that of the drift system: the reduction, by the normal form symmetry of the 
right-hand action of the torus generated by k G so (3), of the left invariant drift 
Hamiltonian tt >— ► tt • £ e + i7T*iV Q 7r on T*SO(3) with canonical symplectic form. 
That the drift Hamiltonian is invariant under this toral action is predicted 
by the theory and is evident by the equality J mi — Jn2- This flavor of the 
theory is subject to the caveat the the frequencies of the linearization of the 
equilibrium are not in 1-1 resonance with the rotation frequencies of the relative 
equilibrium itself; were there to be such a 1 -1 resonance t he preq would acquire 
the structure of a magnetic moment, as in Patrick 199^1 ■ For the case N = 3 



this is not an issue since the reduced spaces are points and so there are no 
reduced frequencies. For N = 4 one verifies by calculating the linearization that 
if the reduced frequencies are ±w re( j then 

l£e| 2 Q (9 cos 2 a + 4 cos a + 3) sin 2 a 



3 cos 2 a + 2 cos a + 3 

Since the second term does not vanish for < a < tt, the group frequencies ±£ e 
are never equal to the reduced frequencies and the simpler perturbation theory 
suffices. 

If the location of the preq on the sphere is denoted by y, then the equations 
of motion for it (i.e. the equations of motion of the drift system) turn out to be 

d2 y M 2 y , D d y u i , Qfi x 

m a ^ = -m a — -+*Bxv, v = -, B = ^y, (26) 

where a the momentum associated to the normal form symmetry and 

1 



JniR 2 ' 

Thus, the drift system is the same as that of particle of mass m a and charge a 
moving on the sphere under the influence of the magnetic monopole B. Given 
that one is perturbing a relative equilibrium where the central vortex has been 
located at Hk, the initial location of the preq may be taken to be j/(0) = 
Kk, while its initial velocity v(0) and the charge a may be obtained from the 
momentum perturbation Afi by the equations 

m a K z 



By solving Equations (E6J), one finds that the prediction of the theory is that the 
preq will rotate about the perturbed momentum, say A/i, with angular velocity 
JjviA/i. 



Specified momentum perturbations of the relative equilibria p e may be ac- 
complished by moving the central vortex from its original location at i?k to 
(RsiaS, 0, RcosS) while changing the angle a to a + Aa, as follows. The mo- 
mentum of the perturbed configuration is easily verified to be 



JS p 



>-l)I\ 







cos(a + Aa) 



sin 5 



cos 5—1 



so this is A/i = [A/ii, 0, A//3] if 



sin 6 



Am 



cos(a + Aa) — cos a 



(N - 1)F : 



■(A/i 3 + r(cOS<J-l)), 



(27) 
(28) 



and arbitrary directions in A/i may be obtained by rotating this. As these for- 
mulas are to be used to perturb a relative equilibrium it is understood that A/j,i 
and A/X3 are small; just how small depends on the validity of the approximation 
that is the drift system. The momentum of the relative equilibrium, being zero, 
does not itself provide a scale. Certainly, though, the perturbation should not 
significantly affect the geometry of the relative equilibrium, meaning it should 
not displace the vortices by amounts comparable to the diameter 2a. In the 
using of (|27|) and (Eq) one should therefore ensure 



«1, 



Aa 
a 



«1 



(29) 



Using (p7|) and (p3) I have simulated the 4- vortex preq corresponding to a = tt/6 
on the unit sphere for momenta A/x small multiples of the vector [2,0,3], and 
calculated the angular velocities of the preq. The results are on the right of 
Figure (|l|). The drift rates fit well to a cubic with slope | J41 1 at zero. 

An essential aspect of the preq is that their location is ill defined as a concept. 
Suppose for example one assigns the location Rk to a preq corresponding to 
the state po £ P sp at some particular time, and sometime later the state of 
the system is p\ £ P sp . If there is a group element A £ SO (3) such that 
Apo = p\ then the location corresponding topi is RAk, unequivocally. However, 
for N > 4 the the reduced phase space for the vortex system has dimension 
2N — 4 > 4, so the flow on this reduced phase space is usually at least as 
complicated as a toral flow with two incommensurate frequencies. Thus the 
reduced flow may never repeat itself, and there may never be an A such that 
Apo = Pi exactly. For the simulations above the location of the preq was 
decreed to be the average of the locations of its constituent vortices; another 
possibility for example is the location of its central vortex. The preq alters its 
shape on a fast time scale compared with its overall motion, and this gives a 
statistical character to the meaning of its "location" on the sphere. The location 
of the preq becomes ever more exact and its character ever less statistical as the 
perturbation falls to zero. 
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Figure 1: Left: A A-vortex preq with opening angle a = it/ 6 subjected to a large perturbation 
subsequently moves on the sphere in a circular path about the perturbed momentum, which 
points towards you. The sense of rotation as is clockwise, the reverse of the right hand sense 
obtained from the perturbed angular momentum. Shown are the paths of the central vortex and 
one outer vortex; the paths of the other two outer vortices are suppressed to avoid cluttering 
the picture. Upper right: the drift rate vs \Ap,\; the curve shown is a cubic with slope at 
equal to J41. Lower right: the angular velocity of a preq divided by |A/u[, calculated over 200 
consecutive intervals of its path. 



This problem is illustrated in the bottom right of Figure ([l]) , where path of a 
preq has been regularly sampled and plotted are the values of J41 from angular 
velocities obtained from consecutive changes in the angle that the preq makes 
with its initial condition. Were the preq to exactly follow the drift system the 
plot would be constant, but instead one gets a variation about the average value 
of 0.4602. This average is just .07% off the value of J41 = .4599 predicted by 
Equation (|4j). 

The issue of whether it is possible to define the location of a preq may be put 
more deeply as follows. Let the Marsden-Weinstein reduced space for the vortex 
system at its zero momentum level be Pg P with reduced Hamiltonian H sp ; the 
relative eq uilibrium p e c orresponds to the equilibrium, say p e of the reduced 
system. In Patrick |l995 it is shown that there is a symplectomorphism defined 
near the group orbit 50(3) • p e and onto a neighborhood of p e times the zero 
section of T*SO(3) such that the Hamiltonian H sp becomes 



H sp (x, tt) = H sp (x) + -k ■ £ e + r-^TT 2 + h.o.t. 

Thus the reduced degrees of freedom are linked to the overall motion by the 
higher order terms depending on both x and n, and if these terms vanish then 
the location of the preq may be unequivocally assigned by following the fiber of 
T* SO (3) to its zero section. Conversely, the location of the preq is ill defined 
in as much as one is obstructed in removing these higher order terms. 



3 Small Opening Angles 

With opening angles as large as a. — n/6, such as in Figure (fil), the relative 
equilibrium p e does not have the appearance of a localized particle. But one 
may choose a arbitrarily small, and now I discuss some aspects of small a and 
report some numerics for a = 7r/2048 « .088 deg. 

Firstly, in the a — > limit there is an essential difference between the 3- 
vortex and 4- vortex preq: for the 3- vortex preq, from (|22j), 

m a = —8ir + 0(a 2 ), 

while for the 4- vortex preq, from (E4j) , 

Q 

m a = —tto 2 + 0(a ). 

o 

Since the mass of the 3- vortex preq does not fall with a while that of the 4- vortex 
preq does, the 3- vortex preq is very heavy in comparison to the 4- vortex preq for 
small a. Thus, large momentum perturbations are required to move a 3-vortex 
preq. As shall be shown immediately, momentum perturbations must also fall 
with a, or else destroy the relative equilibrium by too grossly perturbing it. 
Thus, for dynamical purposes the 3-vortex preq is infinitely heavy in the a = 
limit. This is confirmed upon passage to the planar system, as will be seen in 
Section (4). 

The restrictions ( |29| ) in the a = limit have a different character as they 
affect A/ii vs. A/13. The effect on A/^i is clear: directly from ( |27| ) and the first 
of © 

I Ami I = |rsin£|<0(a). 

However, by elementary manipulations, (Eq) becomes 

|Aa| = I cos _1 (zcosa) — a\, 



where, temporarily, 



Z = CC . S( 5+^. (30) 



By considering the function z i— > | cos ^zcosa) — a| for small a one verifies 
that J cos -1 [z cos a) — a\ <C a is equivalent to 



2-1 \-Z 

< -^ << 1 or < ~ << 1. 

1 1 cos la 1 

cos a cos a 



In any case this amounts to \z— 1 1 < 0(a 2 ) and since i<a one gets by (pfl) that 
A/i3 <C 0(a 2 ). The point is that for momentum perturbations A/Xi, A/^3 to be 
considered small the first must fall as a while the second must fall as a 2 . Since 
|A/ii| 3> IA/X3I causes preq motion that is rotation about an axis perpendicular 
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Figure 2: Left: A preq incident from the lower left approaches its anti-preq incident from the 
lower right; a repulsive force is evident. Right: the same preq and anti-preq collide head on 
and break up into dipole pairs. The bars at the lower left indicate the diameter of the preqs. 



to the direction of the preq, the effect is that the preqs can have momentum 
0(a) only for motions 0(a)-close to great circle paths. To obtain preq motion 
along smaller circular paths on the sphere one must take A/ii ss A^3 = 0(a 2 ), 
which implies a much smaller 0(a 2 ) momentum. 

Since the maximum reasonable momentum perturbations are 0(a) and the 
4- vortex mass is 0(a 2 ), the maximum 4- vortex angular velocity on the sphere 
is 0(l/a). Thus greater velocities are available to smaller opening angles; the 
velocity available to a 4-vortex preq compared with its diameter is 0(l/a 2 ). 

A most interesting aspect of preqs with small opening angles is that more 
than one of them may be positioned on the sphere, all initially far apart (com- 
pared with their radii) from one another. Since the vortex-vortex interaction 
falls as the the vortex-vortex distance increases, as long as the separate preqs 
remain separated they will not strongly interact with one another, and they will 
maintain their separate identities. However, as they separately move on the 
sphere they may closely approach one another. When two or more preqs closely 
approach the above perturbation theory becomes inapplicable and they undergo 
an interaction, and may move apart from one another largely unchanged, or may 
be partially or completely destroyed. There is no theory of these interactions 
at this time but the main features, obtained by numerical simulation and illus- 
trated in Figures (||) and (||), are as follows: 

• Two preqs interact when their distance falls to lengths comparable with 
their diameters. 

• Two preqs with vortices having the same sign are attracted to one another 
while if the vortices have opposite signs they repel one another. Thus 
identical preqs attract while a preq and its anti-preq (obtained from the 
first by changing the signs of each vortex) repel. 

• The collision of two identical preqs, partially since those preqs attract, 
tends to result in two groups of four vortices representing a state of the 
vortex system far away from the relative equilibria p e . In this sense the 
4-vortex preqs are usually destroyed in same-sense vortex preq collisions. 
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Figure 3: Left: A preq incident from the lower left approaches an identical preq initially at 
rest in the center; and attractive force is evident. Right: a preq moving rightward and slightly 
upward collides with an identical preq moving leftward and slightly downward. The denser 
tracks at the center are the incident preqs. The preq incident from the left remains and 
rebounds back the way it came. The preq incident from the right has been nearly destroyed 
and is a four vortex motion far from the original relative equilibrium. Clearly visible in the 
interaction is the exchange of an outer vortex. 



In such an interaction the two 4- vortex preqs usually exchange one of their 
outer vortices. 

• A very energetic collision of a 4-vortex preq with its anti-preq usually 
results in the destruction of both preqs into vortex dipole pairs. 

• The detailed results of an energetic or non-elastic 4-vortex preq collisions 
are extremely sensitive to initial condition, while the gross aspects of the 
collision (e.g. whether or not the collision is elastic) are relatively more 
robust. 

Many a happy hour may be spent watching the antics of 4-vortex preqs colliding 
of the sphere and the nature of the 4-vortex interaction is the subject of ongoing 
investigations. 



4 Transcription to the Planar System 

As already noted (see Equations (g)), when restricted to sufficiently small re- 
gions of the sphere, the system of point vortices on the sphere reverts to the 
system of point vortices on the plane. This is pertinent for the study of small- 
opening-angle preq interactions, since the interactions occur only as the preq 
closely approach. Moreover, the planar case, being an example with noncompact 
symmetry group, is currently interesting, as the general focus of the Hamilto- 
nian stability literature has been on the compact case, while the noncompact 



situation is known to have distinctive features, as shown by Leonard and Mars 
den [1997]. Consequently, I close this article with a short transcription of the 
above the the planar case; the results should be viewed as the a = limit of 
the spherical system. 

The transcription is easy: through the pull back ([?]) of the sphere to the 
plane, the relative equilibria p e with opening angle a have radii a = Rsina, 
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and by (EG) they rotate at angular velocity 



T / 1 \ N T 

o 2— [iTr — T+cosa )=— — --— ^+0(a 4 ), (31) 

AirR 2 sin 2 a \ N - 1 J N - 1 Ana 2 x h K J 



while, from (|l~3|), 



— r — r 

Ti = \ = — + 0(a 2 

(JV-l)cosa N-l y 



Deleting the higher order terms and renaming a to a, one expects planar relative 
equilibria with radii a, central vortex strength T, outer vortex strength T/(N — 
1), and se(2) generator (8 a ,0), where 

N r 



N - 1 Ana 2 ' 



and it is easily verified that these are indeed relative equilibria of the planar sys- 
tem. Moreover, these relative equilibria occur for the vortex strength parameters 
r„ such that ^2 r„ = 0, parameters for which, by (|10|) , the momentum map J pl 
is equi variant. 

For my purpose, equivariance of the momentum map is important since, in 
the general, equivariance guarantees, for a relative equilibrium p e with momen- 
tum jjL e and generator £ e , the momentum- commutation relation coad{ e [i e = 0, 



and this is crucial for the perturbation theory of Patrick [1995]. On the other 
hand, in absence of equivariance one has the more complicated commutation 
relation 

coad 5e fi e = -S(£e, ') = -«?e S ) 

and the perturbation theory would have to be extended in a fundamental way 
to cover the nonequivariant case. That not being necessary, one calculates the 
nilpotent parts of the linearizations of the planar system at the relative equilibria 
above, to obtain, for the 3-vortex relative equilibria 

1 



4 















2?ra 4 
while for the 4-vortex relative equilibria the nilpotent part turns out to be 




Going to the Lagrangian viewpoint, and comparing with the Lagrangian of a 
2-dimcnsional rigid body moving in the plane, one sees that the 3-vortex preq 
have infinite mass while the 4-vortex preq have inertia coefficient — ^a 2 and 
mass ^ a 2 . 
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